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Abstract 



The feasibility of a perturbation expansion for Green's functions of the 
t — J model directly in terms of X-operators is demonstrated using the Baym- 
Kadanoff functional method. As an application we derive explicit expressions 
for the kernel G of the linearized equation for the superconducting order pa- 
rameter in leading order of a 1/N expansion. The linearized equation is solved 
numerically on a square lattice taking instantaneous and retarded contribu- 
tions into account. 

Classifying the order parameter according to irreducible representations 
Ti,i = 1, ...5, of the point group C^ v of the square lattice and according to 
even or odd parity in frequency we find that a resonably strong instability 
occurs only for even frequency pairing with d-wavelike symmetry. The 
corresponding transition temperature T c is ~ 0.01 \t\ where t is the nearest- 
neighbor hopping integral. The underlying effective interaction consists of an 
attractive, instantaneous term and a retarded term due to charge and spin 
fluctuations. The latter is weakly attractive at low frequencies below ~ J/2, 
strongly repulsive up to ~ |t| and attractive towards even higher energies. 
T c increases with decreasing doping 5 until a d-wavelike bond-order wave 
instability is encountered near optimal doping at 5 bo ~ 0.14 for J = 0.3. 
T c is essentially linear in J and rather insensitive to an additional second- 
nearest neighbor hopping integral t'. A rather striking property of T c is that 
it is hardly affected by the soft mode associated with the bond-order wave 
instability or by the Van Hove singularity in the case with second-nearest 
neighbor hopping. This unique feature reflects the fact that the solution of 
the gap equation involves momenta far away from the Fermi surface (due to 
the instantaneous term) and many frequencies (due to the retarded term) 
so that singular properties in momentum or frequency are averaged out very 
effectively. 

PACS numbers: 74.20-z, 74.20Mn 
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I. INTRODUCTION 



electronic orbital 
the flavor index 



Many studies of the t — J model suggest that this model is able to describe a large body 
of the low-energy physics of real high-T c superconductorsBi. This is true for many normal 
state properties of high-T c oxides where accurate numerical predictions of the t — J model 
are available for the comparison with experiment. Whether the phenomenon of high-T c 
superconductivity itself can be explained within this model is presently not so clear. There 
are several calculations yielding instabilities of the normal state with respect to d-wave 
superconductivity^ and also reasonably high values for the transition temperature tJHI. 
These calculations, however, use often somewhat uncontrolled assumptions, making definite 
conclusions difficult. There is also the viewS that the large observed values for T c are not 
directly related to large mean field T c 's in isolated C11O2 planes described by the t — J 
model. Instead it is argued that pair tunneling between planes enhances strongly weak, 
plane-related superconducting instabilities producing in this way the phenomen of high-T c 
superconductivity. 

In this paper we present a new attempt to calculate T c for the t — J model in a well 
controlled way. Similar as in Refs.i'§ we do not assume the validity of Migdal's theorem or 
approximate the self-energy by the lowest skeleton graphs. Instead we assume that 1/N can 
be considered as a small parameter where N is the number of electronic degrees of freedom 
per site. N consists of two spin directions times N/2 copies of the loca 
counted by a flavor index. Similar like in many slave boson calculation 
is introduced in a somewhat artificial way just to make N a large integer. The £77(2) spin 
symmetry of the original model is thus enlarged to the symplectic group Sp(N/2). The 
original constraint of having no double occupancies of sites is modified to the condition 
that at most N/2 electrons can occupy the N states at each site. Compared to Refs.ifl our 
treatment exhibits two novel features. First, the constraint is implemented in a different and 
more rigorous way yielding differences in the equation for the superconducting gap already in 
the leading order 0(1/N). Secondly, we have solved the linearized gap equation numerically 
obtaining also values for T c . Our conclusions about the occurrence of superconductivity in 
the t — J model are thus no longer based only on Fermi surface averaged, static coupling 
strengths as in Refs.i'§00. 

Regarding constraints their implementation in the X-operator approach is trivial. 
Mathematically, the constraint means that the sum over diagonal elements of X-operators 
at a given site has to be equal to N/2. This sum commutes with every X-operator and 
thus is a multiple of the identity operator in any irreducible representation of X-operators. 
Enforcing the constraint means therefore just to select the correct subspace of Hilberts space 
where the eigenvalue of this sum is equal to N/2. In slave boson theory the X-operators are 
represented by product of slave operators and the Hilbert space is enlarged to the Hilbert 
space of slave particles. The constraint means now that at each site the number operator of 
slave particles has to be equal N/2. This condition can obviously not hold as an operator 
identity, i.e., cannot be enforced at any position of an expectation value of slave operators. 
However, it should be enforced at all positions which originally separated X-operators. Such 
an enforcement for a fixed N seems to be in conflict with the Bose condensation of the 
bosonic slave particles as well as the independence of fermion and bosonic slaves in the limit 
N —>■ 00. Using the Dirac method to enforce constraints on the operator level it has indeed 
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been shownEl that the commutator relations for the slave particles have to be changed and 
that, for instance, the fermionc and the bosonic slave operators no longer commute with 
each other. In view of this open problems we use in sections II and III a perturbation 
expansion directly in X-operators following the work of Ref.0 This procedure will give us 
also the opportunity to compare the 0(1/N) expression for the gap equation with that of 
the slave boson approach to see whether the 1/N expansions are really the same in the two 
cases. The final analytic results for the gap equation have already been presented in Ref.i, 
however, without derivations. These derivations can be found in sections II and III. 

The second new feature of our work deals with the solution of the gap equation. Pre- 
vious work using 1/N expansions concluded from Fermi surface averaged, static coupling 
constants on the occurrence of superconductivity. Our previous worko0 has shown that 
superconducting instabilities occur for any symmetry channel, doping and both for J = 
and J ^ 0. The point is more whether the corresponding T c is very small or large enough 
to be relevant. T c here is understood as the transition temperature within mean-field the- 
ory. This means that the lowering of T c due to fluctuations in the superconducting order 
parameter is not taken into account. This assumption is justified if one compares with real, 
three-dimensional superconductors but does not apply, of course, to strictly two-dimensional 
models where fluctuations drive the order parameter to zero at any finite temperature. In 
order to determine T c one has to solve the gap equation which has several non-BCS features: 
The kernel of the linearized gap equation consists of an instantaneous and a retarded term 
and both are characterized by different cutoffs. Moreover, the presence of the instantaneous 
term does not allow to restrict the momenta to the Fermi surface. Consequently, we solved 
the gap equation by numerical means. Our method is also suitable to investigate the ques- 
tion of odd frequency pairing|ll in our model. In this case the instantaneous term drops out 
but the full frequency dependence of the kernel must be kept in addition to the momentum 
dependence along the Fermi line. Results for odd frequency pairing will also be presented 
in section IV together with the conclusions. 



II. MODEL AND GENERAL EQUATIONS FOR THE ELECTRON GREEN'S 

FUNCTION 

The Hamiltonian of the t — J model can be written as 

H = E %xfx° p + v ^-xfxf- v ^- r xfxf. (i) 

AT 1 3 Z-^ AfJ * 3 AAT 1-3 V ' 

ij 1 ij ^ ij 

p= 1...N p,q— l...iV p,q — l...N 

Let us consider first the case N = 2. Assuming one orbital per site and excluding double 
occupancies of sites there are three states \ p > at each atom i. p = denotes the empty 
state, and p = 1, 2 singly occupied states with spin up and down. The Hubbard operators 
Xf q can be represented as projection operators Xf q = \ p . >< q A and obey the following 
commutator and anticommutator relations 

[xp, x; s ] ± = %(<v*f ± Ssrxn, (2) 

and the completeness relation 
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p=0 



(3) 



The upper (lower) sign in Eq.(2) holds for bosonlike or mixed (fermionlike) Hubbard opera- 
tors defined by p, q > or p = q = (p = 0, q > or p > 0, q = 0). For i — j both the upper 
and lower signs hold in each case. The first term in Eq.(l) describes the hopping of particles 
between the sites i and j with matrix elements tij. The second term in Eq.(l) denotes the 
Heisenberg interaction between the spin densities at site % and j with the exchange constants 
Jij. The third term in Eq.(l) represents the charge-charge interaction of the t — J model. 
In the following we consider only between nearest neighbors (J^ = J) and between 
nearest (ty = t) and next nearest (t^ = if) neighbors. We also use always \t\ as energy unit. 

The Hamiltonian in Eq.(l) is a generalization from N = 2 to an (even) arbitrary integer 
N. The orbital index p consists now of the spin and a flavor index, the latter enumerating 
N/2 identical orbitals at a site. The symmetry group of H is the symplectic group Sp(N/2). 
For N > 2 the operators X can no longer be written as projection operators in some 
basis. Instead, fermionlike (bosonlike or mixed ones) operators are assumed to satisfy the 
anticommutator (commutator) relations of Eq.(2). Only some of the diagonal operators can 
assumed to retain their projection properties, namely, (Xf p ) 2 = Xf p for p > 0. Most other 
relations characteristic of projection operators such as X^X® 1 = X} 1 are lost for N > 2. 
The completeness relation Eq.(3) is replaced by the constraint 

N 

Qi = Y. x f P = N / 2 - ( 4 ) 

p=0 

By explicit construction of the Hilbert space and the action of the X's on its vectors one 
can show@ that the above properties, together with H, specify completely the problem. In 
particular, Xf° is a non-negative operator. As a result, Eq.(4) means that at most N/2 
particles can occupy the iV available states at a site. In this way one may expect that 
expectation values of observables approach smoothly the physical case N = 2 from large 
N's yielding a basis for 1/N expansions. We also would like to point out that slave boson 
treatments of H have many features in common with our approach. However, the Hilbert 
spaces are different in the two cases as well as the enforcement of the constraint Eq.(4). 
Qi commutes with all Hubbard operators and thus is proportional to the identiy in any 
irreducible representation of the X-operators. Enforcement of the constraint means in our 
case just the selection of the correct subspace of the Hilbert space where Eq.(4) is satisfied 
as an operator identity. 

Following the Baym-Kadanoff procedureEUli we define a non-equilibrium Matsubara 
Green's function of two fermionic operators X(l) and X(2) by 

G(12) = - < TSX{1)X(2) > / < S >, (5) 

s = Te fdlX(l)K(l) (6) 

The number 1 stands for an internal pair index p, q as well as a site and a (imaginary) time 
index i,r. J dl means J2 P , q ,i So where (3 is the inverse temperature. In Eqs.(5) and (6) 
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T is the time ordering operator and K an external source term which is assumed to couple 
only to bosonic X-operators. We also introduce the non-equilibirum expectation value of 
bosonic X-operators by 

L(l) =< TSX{1) > I < S > . (7) 

Using the Hamiltonian Eq.(l) the Heisenberg equation of motion for fermionic operators 
becomes 



_d_ 

with 



X(l) = J c/2d3t(123)X(2)X(3), (8) 



£(123) = 5(r 2 - n)5(r 3 - r x ) - n * 3 1112 + N * 112 ^ — S qi0 S q30 (l - S P30 )5 P20 5 q2qi 8 P3Pl + 5 P2P1 8 q2P3 ) 

— ^pi0^p 3 o(l — ^g 3 o) (^293^291 + dpiPl ^320^391 ) • (9) 

Eq.(9) agrees with Eq.(5) of Ref.El if the different sign convention for the hopping term is 
taken into account. Using Eq.(8) and rewriting higher-order correlation functions in terms 
of functional derivatives with respect to K the equation of motion for G can be written as 

J d2(Go X (12) - E'(12))G(21') = Q(ll'), (10) 

E'(12) = -J d3£(132)L(3) + J rf3rf4rf5t(134)G(45)r(52, 3), (11) 

T(12, 3) = 5G- 1 (12) /5K{3). (12) 
Go is the unperturbed Green's function 

^(12) = -6(1 - 2) A _ 5(1 - 2)(K 00 (l)5 qiq2 - K™{1)). (13) 

Q is given by 

Q(ll') = 5(1 - l')(L™'(l)5 qp , + LP\l)5 q/p ), (14) 

where the index 1 = (pq,i,r) has been split into 1 = (pq,l) with 1 = (ir). L can be 
expressed by G so the above system of equations for G, the self-energy S , and the vertex V 
is closed. 

For our purposes it is more convenient to use a normalized Green's function g which 
obeys a Dyson equation with the usual 5-function on its right-hand side. Writing 

0(110 = / d2G(12)Z- 1 (2Y), (15) 

and requesting that Z satisfies the following equation 
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Z(W) = Q(ll') - J d 3dAd5t(13A)g(A5) 6 -^^-, (16) 

Eqs(10)-(12) assume the form 

J d2(G \l2) - £(12))< ? (2l') = 5(1 - 1'), (17) 

£(12) = - J d3*(132)L(3) + J d3d4d5t(134)#(45)7(52, 3), (18) 

7(12,3) =5g-\l2)/5K(3). (19) 
From Eqs.(17) and (19) follows, moreover, the equation for the vertex 

7(11', 3) = a(W; 3) + J dAd5Q(ll', 45) 7 (45, 3), (20) 

with 

5G \\Y) t^ SZjlV) 5L(A) 

a{n ' 3) = -JkIsT ~ J dA ^W • Ik®' ( } 

0(11', 45) = / ^7^^^(64)^(57). (22) 

In Eq.(22) we have replaced the functional derivative of g with respect to K by the vertex 
7 using Dyson's equation which also leads to a sign change. The above equations are exact 
and hold both in the normal and the superconducting state. In the normal state the two 
internal index pairs in G(12) must obey either pi = q2 = or p2 = q± = 0. From Dyson's 
equation follows then that the same holds for Q(12) whereas the possible indices in G (12), 
£(12), and #(12) obey either pi = p 2 = or qi = q 2 = 0. In the superconducting state there 
are no such restrictions; the only general requirement is that all these indices are associated 
with fermionic operators. 

The self-energy is a functional of g and L, i.e., 

£ = £[<?,£]. (23) 

In order to derive the linearized equation for the anomalous self-energy we split £ into the 
normal part £jv and a small anomalous part £ an . Expanding Eq.(23) up to linear terms we 
obtain 

E » (11,) = / d2d3 (fj|y)^( 23 ) + Mw)/ L < 2 >' < 24 > 

where 5g and 5L are of first order in the anomalous part and the subscript iV means that the 
functional derivatives are to be taken in the normal state. Now 5L has no linear contribution 
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due to gauge invariance and thus drops out in Eq.(24). Calculating 5g from Dyson's equation 
Eq.(17) in linear approximation we obtain 

E an (ll') = J d2d3e(ll',23)I! an (23), (25) 

where is given by Eq.(22) with all quantities taken in the normal state. Eq.(25) represents 
the linearized equation for a general order parameter for superconductivity. In order to have 
superconductivity the homogenous equation Eq.(25) must have a nonvanishing solution for 
S an . The highest temperature where this occurs is the transition temperature T c . 

The basic quantitity to be calculated is according to Eq.(25) the kernel 6. In order to 
find a functional equation for 9 we get from Eq.(18) 

= / d2t(127) 7 (81', 2) + J d 2d3d6d7d8t(123)g(36) 6 -^^, (26) 

and from Eq.(20) 

57(61', 2) _ <Ja(61',2) 



5g(78) ^ 5g(78) 

/ d9dl0d9'dl&e(61', QlO') ^ 9 ' 10 '' 2 ) + / d9 / dl0 ^ e( r 61 / 9 ' 1 ° ,) 7(9 , 10 / , 2). (27) 
J dg{78) J dg{78) 

Solving Eq.(27) for S^/Sg, inserting the result into Eq.(26) and inserting Eq.(26) into Eq.(22) 
yields an exact functional equation for 0. 



III. 1/N EXPANSION FOR THE KERNEL 6 

The equation for G obtained in the previous section is too difficult to be solved directly. 
On the other hand it is exact and holds for any iV so it may serve as a starting point for 
approximate treatments. In the following we assume that 1/N may be used as a small 
parameter and calculate 6 up to order 1/N. This will allow us to obtain the gap equation 
Eq.(25) in leading order of the 1/N expansion. 

The iV-dependence of an equilibrium quantity is determined by the number of coupling 
constants it contains and the number of free internal summations. In equilibrium, i.e., 
without the source term K, we have, for instance, 

9(°f °| 2 ) = <W(l-2), (28) 
S(°f°| 2 ) = ^ 2 S(l-2), (29) 

l(°f °| 2 , ™ 3 ) = -^p 3 ?3<W7c(i2, 3) + 5 qiP3 5 q2q3 -f s (12, 3), (30) 

and similar expressions for a and 6. In Eqs.(28)-(30) all the indices qi,q2, ■■■ are assumed 
to be larger than zero. The equilibrium Green's function g, the self-energy S, the charge 
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vertex 7 C and the spin vertex 7 S are free of internal indices. In the following we need only 
the leading 0(1) contributions for these quantities. From Eqs.(13), (17) and (19) follows 
that the spin vertex becomes then simply equal to 5(1 — 2)5(1 — 3). Using the constraint 
Eq.(4) one also recognizes that the charge vertex 7 C defined in Eq.(30) is equal to the element 
<?i = <?2 > 0,p 3 = g 3 = of the general vertex 7 on the left-hand side of Eq.(30) which is the 
motivation to use prefactors in Eq.(30) in defining 7 C . 

Decomposing the labels 1,1', etc. into their internal and external parts the anomalous 
self-energy or order parameter has the form S an ( p i° ), or in more detail, E^^ 1 ™ 1000 ^™ 1 )- 
In the following we are interested in order parameters which are structureless in the flavor 
indices, i.e., we assume always m 1 = m[. With respect to spin indices we consider either 
singlet or triplet pairing. Both cases are included if we put U\ = —cr[. Correspondingly we 
will write for the internal indices pi,pi thus indicating the special relationship between these 
two indices. Exposing explicitly the internal indices the gap equation Eq.(25) becomes 

sw*! ? ) = E / ^ , f f )E.„ ( f f), pi) 

where O is to be taken in the normal state. Eq.(31) also makes visible which indices for 
O are actually needed in the gap equation. A closer examination shows that the first two 
contributions on the right-hand side of Eq.(27) are of higher orders in 1/N than the third 
term in this equation for the needed combination of indices.. Dropping these terms and 
inserting Eqs.(26) and (27) into Eq.(22) yields the following equation for O valid up to 
0(1/N) and for the above combinations of internal indices: 

0(11', 910) = J d2d7tZ8t(127)7(81', 2)0(7%(1O8)+ 

r 50(61' 910') 

/ d2d3d6d7d8d9'dl0't(123)g(36) — \, / . ; 7 (9'10', 2W79W108). (32) 
J 5g(7S) 

The first term on the right-hand side of Eq.(32) could have been also obtained directly from 
Eq.(18). Considering only anomalous contributions to S the first term on the right-hand 
side of Eq.(18) cannot contribute, wheras the second one can contribute in two ways: a) 
via an anomalous g and a normal vertex 7, b) via a normal g and an anomalous vertex 7. 
Case a) corresponds to the usual situation of Eliashberg theory and also of the slave boson 
treatment of superconductivity in the t — J model in 0(1/N): S ara consists then of a Fock 
diagram containing an anomalous Green's function and an effective interaction taken in the 
normal state. Linearizing g in S ara case a) immediately yields the first term on the right- 
hand side of Eq.(32). Its evaluation is straightforward and shows that it is of 0(1/N) for 
the relevant combination of internal indices. This indicates that T c — > for N — > 00 which 
agrees with the fact that at N = 00 the system consists of renormalized, but non-interacting 
fermions. 

The leading contribution of the second term on the right-hand side of Eq.(32) is obtained 
by taking the spin-flip part in the hopping matrix element t and the spin vertex 7 S in 0(1) 
for the vertex 7. As a result only one sum over internal indices survives. The second term in 
Eq.(32) is of 0(1 /N) if the two functions 
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# (1) (61',2|78) 



E 

P3 



50 ( 



P3O Opi 
6 V 



22/ 



Qpi 
8 



(33) 



m / P3U Upi pjjU piU \ 

^(61-,2|78) = E ■ (34) 

are of 0(1). Equations for i — 1, 2 can be obtained from Eq.(32) by taking appropriate 
derivatives. An examination of the various terms shows that terms with a second functional 
derivative of are smaller by a factor 1/N compared to the leading ones and thus can be 
neglected. One then finds that the functions have the form 

9 (1) (61',2|78) =6(2- 8)^(61', 2|7), (35) 



2 (2) (6T,2|78) =5(2-7)(7 (1) (61',2|8). (36) 

The reduced functions grW(61', 2|8) (the use of the same name for corresponding functions 
with different number of arguments should not cause any confusion) satisfy the integral 
equation 

0®(l,r,2|7) = fcW(l,r,2|7) -|d3c/4d6^t(l3%(64)^)(6r,3|7)^(23), (37) 

with 

/i (1) (ll', 2|7) = M(Il , 7)^(21'), (38) 



/i (2) (ll',2|7) = - / d3c/4iVt(13 , 4)7 c (7l', 4)^(23). (39) 

Here t (123) denotes the external part of £(123), Le., the first three factors outside of the 
bracket in Eq.(9). Inserting the explicit expression for t into Eq.(37) one finds that Eq.(37) 
represents an integral equation for with a kernel consisting of not more than 6 separable 
contributions. The solution of Eq.(37) thus reduces to the inversion of at most 6x6 matrices 
similar as in the case of the charge vertex j c dicussed in Ref.0 Inserting the resulting 
expressions back into Eq.(32) and performing Fourier transforms one obtains after some 
elementary, but tedious algebra the following results. The linearized gap equation Eq.(25) 
becomes 

=-(*) = -]^i:e(*,V)^L_ J E.(V). (40) 

iV c is the number of cells and k the supervector k = (n, k), where n denotes a fermionic 
Matsubara frequency u n = (2n+l)7rT with T being the temperature. e(k) is the one-particle 
energy with momentum k in the limit iV —>■ 00 given by the implicit equation 

c(k) = 6 -t(k) - J(k) ■ 1- £ cos(p x )f(e(p - fi)) (41) 
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where 5 is the doping, t(k) and J(k) the Fourier transforms of the coupling constants tij 
and Jij, respectively, fi is the chemical potential and / the Fermi function. 
The kernel G in Eq.(40) consists of four different terms 

S(k, k') = Q (1 \k, k') + (2) (A;, k') + (3) (A;, k') + 6 (4) (A;, k'). (42) 

The first two terms are given by 

QW(k, k') = Tt(k') t J(k - k') - t(k') - J(k - k'), (43) 

0( 2 )(A;, k') = (i(k') + J(k - k , ))( 7c (Jfe, k' - k) + 1). (44) 
Here, 7 C is the charge vertex in 0(1) determined by 

lc (k,q) = -l+ ]T F a (fc)(l + X (g))^(?), (45) 

a,/3=l 

with the susceptibility matrix 

wW .± gWfW g^gg (46) 

g is the supervector g = (q, iz/ n ) where z/„ denotes the bosonic Matsubara frequency 2nnT. 
The two vectors E and F are given by 

E a (k, q) = (1, i(k + q) + J(q), cosk x , sink x , cosk y , sink y ), (47) 

F^(k) = (£(k), 1, Jcosk x , Jsink x , Jcosk y , Jsink y ). (48) 
The explicit expressions for the third and fourth contributions to are 

0(3)(A;, fc') = - E ^(-k)x2r(A; - fc')7(fc, fc' - *0 + E E s (-k)x rs (k - k')V r (k - k'), (49) 

r=l r,s=l 

Q^(k, fc') = ± E (^(-k)Xrs(A; + fc'))(^r(-k') + ^(A; + fc')), (50) 

r,s=l 

with the vectors V and J¥ 

K(& - hf) = E(l + X(* " k'))^X2, s {k ~ k')lc(k, k'-k), (51) 

s=l 

w r (k + k') = - E (i + x(* + fc'))- 1 ^^-^)^^ + *;')■ (52) 

s,t=l 
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The vectors E and F have five components and are obtained from the six-component vectors 
E and F by omitting the second component. In an analogues way, \2,s is obtained from the 
second row of the matrix \ by omitting the second element. Similarly, the 5x5 matrix \ is 
obtained from the 6x6 matrix x by dropping the second row and second column. The upper 
(lower) signs in the contibutions 6^ refer to singlet (triplet) pairings. 

The kernel Q(k, k') is invariant under the transformations of the point group C iv of the 
square lattice. As a result T, an {k) transforms corresponding to one of the five irreducible 
representations Tj of C^v, and it can be chosen to be either even or odd in Matsubara 
frequencies. Ti corresponds to singlet s- or extended s-wave, T 3 to singlet d-wave- and T 5 
to triplet p-wave pairing for even frequency pairing. The linearized gap equation Eq.(40) 
splits completely into its irreducible parts. This also means that the momenta k, k' can 
be restricted to 1/8 of the Brillouin zone which greatly simplifies the numerical solution of 
Eq.(40). In the following we will deal with even frequency pairing unless the opposite is 
explicitly stated and drop for simplicity the index which differentiates between even and 
odd frequency-pairing. The irreducible kernels and order parameters will thus be denoted 
by Qi(k,k') and T, anji (k), respectivley. 



IV. DISCUSSION OF THE GAP EQUATION AND NUMERICAL RESULTS 

According to Eq.(42) the total kernel is the sum of four contributions. The first one, 
0W, represents the instantaneous part of 6. From Eq.(40) follows that positive values for 
@(!) correspond to repulsion, negative values to attraction between electrons. Keeping only 
and decomposing it into its irreducible symmetry components the gap equation Eq.(40) 
becomes for the representation i=1...5 of the point group C iv of a square lattice 

Si,«(k) = -^Ee, (1) (k-k')x(k')^(k'), (53) 

with 

_ tanh{{e{k') - fi)/T) 
X[k) ~ 2(e(k') - /i) ' (5 j 

e?> = ei 1} = e?> = o, and 

e?(k,k') = 8\t\ Vl (k') - SJ^k^k'), (55) 



Q { 3 l) (k,k') = -8J m (k)r ]3 (k>). (56) 

In Eqs.(55)and (56) we assumed nearest-neighbor hopping with tij = —\t\ and used the basis 
functions ^i(k) = (cos(k x ) + cos(k y ))/2, r]s(k) = (cos(k x ) — cos(k y ))/2. Inserting the kernels 
Eqs.(55) and (56) into the gap equation yields the following conditions for a finite T c : 

/ \ i i ah\ \t \J 
s — wave pairing (E an ~ constant) : 1 + a\t\ + — = 0, (57) 
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extended s — wave pairing (S ara ~ ^i(k)) : 1 — b\ + ^ ^ L = 0, (58) 



d — wave pairing (S an ~ ^(k)) : 1 — b 3 J = 0. (59) 

Here we used the following abbreviations 

8 



NN C ^ 



Ex(k>(k'), (60) 



^ = ^Ex(k')^(k'), (61) 



1 c 



for / = 1,3. A numerical evaluation shows that in the interesting parameter region < J < 
0.3,0 < 5 < 0.8 (5 is the doping) and at low temperaturs a,bi,b 3 are positive and b\ < 3. 
This means that there is never an instability with respect to constant s-wave pairing, bi is an 
increasing, 6 3 a strongly decreasing function with increasing S and these two functions cross 
around 5 ~ 0.6. Thus d-wave pairing is stable in any case for 5 < 0.5. For 0.5 < 5 < 0.8 
a\t\ is much larger than 1 so that the third term in Eq.(58) cancels the second term to 
a large extent making extended s-wave pairing also in this region unfavourable. Thus we 
find that the instantaneous contribution to the kernel strongly favours d-wave pairing in 
the interesting parameter regime. The reason for this is to a large extent a band structure 
effect: 63 is much larger than b\ at small dopings because of large contributions around the 
X-point. At larger doping the competing extended s-wave pairing is strongly suppressed by 
its coupling to the hopping term, which ultimately is caused by the constraint. 

@( 2 \ 0®, and G^ 4 -* are retarded contributions to 0. is mainly determined by collec- 
tive charge fluctuations due to the poles of 7 C . and G^ 4 -* originate from the anomalous 
part of the vertex and involve both charge and spin fluctuations. The latter dominate in 

e(4) 

in agreement with the sign change between its singlet and triplet contribution. Quan- 
titatively, G^ 4 -* is by far the largest of the retarded contributions. Our Q does not contain 
terms which are related to the magnons of the undoped case. Such contributions are of 
higher order in 1/N than those considered above and are thus neglected. 

Our expression for O is different from that of the slave boson approach!. S an in the 
slave boson approach essentially consists of a spinon Fock term with an effective interaction 
taken in the normal state times an anomalous spinon Green's function. Such an contribution 
clearly corresponds to the case a) discussed after Eq.(32), i.e., to the first term on the right- 
hand side of Eq.(32), or, equivalently. to the sum of 6^ and O*- 2 -*. Our contributions O*- 3 - 1 
and ©( 4 ) have no analogue in the slave boson approach though they are clearly also of 0(1/N). 
The differences between the two approaches can be made more explicitly by considering two 
limiting cases. For J->09 reduces to Eq. (9) of Ref .0 The underlying interactions are still 
spin-dependent as can be seen from the singlet versus triplet case. Such a dependence on 
spin is not present in the corresponding slave boson expression^. Furthermore, the high- 
frequency limit of O, i.e., should be identical to that of the slave boson expression if the 
two approaches are equivalent. The slave boson result for O^ differs, however, from Eq.(43): 
the argument in one of the two hopping terms is replaced by k and there is an additional 
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term proportional to 1/5. Eq. (5) agrees, however, with previous X-operator results based on 
the diagrammatic methods for X-operators, for instance, Eq.(35) of Ref.0. The expression 
for Q« of Refl contains only one of the two hopping terms t(k') of Eq.(5) which, we think, 
is incorrect. The general question then arises why is it possible to obtain different 0(1/N) 
expressions for the same quantitiy in the two approaches using the same Hamiltonian. The 
differences arise because different Hilbert spaces are used in the two cases. This difference 
is already indicated in the definition of the order parameters for superconductivity. In the 
X-operator approach all vectors of the Hilbert space are eigenvectors of the constraints Qi 
with eigenvalue N/2. Expectation values of operators are only nonzero if, in the slave boson 
language, the number of created slave particles is equal to the number of annihilated slave 
particles in these operators. The superconducting order parameter of slave boson theory, 
on the other hand, consists of the expectation values < b > and < c'c* >. Clearly, these 
expectation values have no analogues in the X-operator approach. 



A. Eigenvalues of in the static limit 

In the weak-coupling case it is often assumed that the solution of the gap equation 
Eq.(40) is of BCS-type for each symmetry Tj, i.e., 

T ci = l.Uu c e l ' x \ (62) 

uj c is a suitable cut-off and Aj the smallest eigenvalue of a matrix which is essentially equal to 
the static limit of the kernel G and given by Eq.(lO) of RefR This matrix includes also the 
prefactor 1/N on the right-hand side of Eq.(40) and we put from now on N always equal to 2. 
Actually we will find that Eq.(62) is a rather poor approximation for T c because G contains 
instantaneous and retarded contributions which have different cut-offs. Nevertheless the 
study of the static kernel G and its lowest eigenvalues has traditionally played a great role 
in discussing superconductivity in t-J models. 

Fig. la) shows the terms and G^ 2 -* and Fig. lb) the terms 0( 3 ) ; @( 4 ) and for singlet 
pairing at zero frequencies for a fixed first momentum k = (2.465,0.309) as a function of the 
second argument k'. k' moves counterclockwise around the Fermi line passing through the 
points X, Y(X, Y) along the positive (negative) x- and y-axis, respectively. Note that the 
letters X, Y denote not the k-points (?r, 0), (0, 7r) but the points on the Fermi line between 
the T-point and the points (vr,0), (0,7r), respectively. The doping is 5 = 0.17 and J = 0.3. 
We used 22 k-points along 1/8 of the Fermi line and a net of 300x300 k-points in the Brillouin 
zone. Positive values of G mean repulsion, negative ones attraction between electrons in the 
s-wave channel. and originate from the normal, G^ and Q^ 4 -* from the anomalous 
vertex. is due to instantaneous, the other terms due to retarded interactions. The 

Figures show that G^ 4 -* and 0^ are by far the largest contributions. Both are dominated by 
the d-wave component and add up in a coherent way. has different signs for singlet and 
triplet pairings. This as well as the explicit calculation shows that important contributions 
to G*- 4 -* come from spin fluctuations within the partially filled band. There are no spin 
contributions related to the Heisenberg term and the magnon spectrum at zero doping 
because these terms are at least one order in 1/N smaller than the considered ones. 

(3) 

is due to charge (and spin) excitations and it is very small. Finally, Q*- 2 -* is due to charge 
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fluctuations and consists of an attractive s-wave component (which, however, is cancelled 
by a repulsive s-wave term in 6^) and a d -wave part which cancels partially that of the 
term. Figs.l demonstrate the importance of anomalous vertex contributions to 0, the 
d-wave character of the leading contributions, and the competition between instantaneous 
and retarded interactions. 

Fig. 2 shows the lowest eigenvalues Aj for each of the five representations IY In the 
calculation we used 5 k-points along 1/8 of the Fermi line and a net of 300x300 k-points 
in the Brillouin zone. The eigenvalues Aj decrease with decreasing doping 5 and diverge 
at 5 ~ 5 bo ~ 0.13. At this doping value one of the six eigenvalues of the 6x6 matrix 
1 + x in Eq.(45) goes through zero . As a result there is a soft mode which freezes into a 
static incommensurate bond-order wave of d-wave symmetry for 8 < 8bo- This instability 
causes a divergence in some of the contributions to O at 5 bo and, as a precursor, the large 
negative values seen in Fig. 2 for 8 < 0.15. The solid line in Fig. 2 exhibits the lowest 
eigenvalue A 3 associated with d-wave pairing. Its absolute magnitude is much larger than 
all the other eigenvalues and this is true for all dopings. Inserting A, into Eq.(40) it is clear 
that superconducting instabilities for non-d-wave symmetries are more of academic interest 
because the corresponding T' c s would be extremely small. One concludes that the normal 
state is generally unstable against superconductivity in every symmetry channel but that 
only the d-wave (r 3 ) instability is strong enough to account for high transition temperatures. 

B. Frequency dependence of 6 

For either even or odd frequency pairing the frequency dependence of can be rewritten 
in terms of one frequency argument u n which is equal to the difference of the two original 
frequency variables. In order to illustrate the frequency dependence of the retarded part of 
O, rei , in the case of d-wave scattering, we average the momenta in ret over the Fermi line 
in the irreducible Brillouin zone using the eigenvector of the lowest eigenvalue A3. Finally 
we perform the analytic continuation iuj n — > uj + irj. 

Fig. 3 shows the negative imaginary part of ©^(cu + irj) for rj = 0.005. It is the analogue 
of the familiar function a 2 F(u) of Eliashberg theory for phonon-induced s-wave supercon- 
ductivity. In our case this function is no longer positive-definite. This is a result of the 
projection of the total kernel on the T 3 symmetry. This projection involves a sum of k' 
over the 8 symmetry-related points along the Fermi line with factors determined by the 
representation i. In the case of d-wave pairing these factors are positive for small and large 
momentum transfers (of the order or a reciprocal lattice vector) and negative for interme- 
diate momentum transfers of about half of a reciprocal lattice vector. For collective density 
fluctuations described by 0^ 2 ^ this means that high-frequency contributions appear with 
a negative and low-frequency contributions with a positive sign in the d-wave kernel. A 
similar behavior is found for the other contributions. As a result -ImO™' is positive for 
uj < uj 2 ~ \t\ and negative for 00 > 002- Using a momentum average with a constant weigth 
of Q ret corresponding to constant s-wave pairing would yield a curve which is similar to 
that in Fig. 3 for uj < uo 2 but opposite in sign for uo > uo 2 . —ImQ r z et describes the spectral 
distribution and spectral weight of spin and charge excitations involved in d-wave scattering. 
This function extends over a wide frequency region of about 2\t\. Its high-frequency part 
is typical for collective charge fluctuations. It also contains substantial spectral weight at 
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lower frequencies exhibiting a rather linear increase in frequency at the low-frequency end 
of the spectrum. Fig. 4 shows the corresponding real part of 0g et . It is weakly negative at 
low frequencies up to about uj 1 ~ J/2, changes then from negative to positive values until 
about oo ~ oo 2 . The retarded interaction between electrons is thus attractive for < u < U\ 
and strongly repulsive for uj\ < uj < uj 2 . Somewhat above uj 2 the real part of <S>l et jumps 
to large negative values and approaches zero from below in the high-frequency limit. The 
inset of Fig. 4 shows the real part of 3 at small frequencies with a larger resolution using 
r] = 0.002. 

C. Transition temperatures 

Figs, la) and b) show that the instantaneous and the retarded contributions to are 
of similar magnitude and often compete with each other. In the calculation of T c they are 
associated with quite different cut-offs. The instantaneous part is characterized by a cut-off 
determined by the width of the effective band whereas the cut-off relevant for the retarded 
part is set by the frequency range where attraction dominates characterized roughly by J. 
In view of these complications we developed a method to solve Eq. (40) directly, avoiding the 
use of pseudopotentials. The only simplification we use is to put the momenta of the retarded 
kernel on the Fermi line. The validity of this approximation has been checked numerically 
and holds very well in our case. Using the fact that the instantaneous kernel consists only of a 
few separable contributions Eq.(40) can be reduced to a linear matrix problem. The number 
of rows and columns are given by the number of considered Matsubara frequencies times 
the number of k-points on the Fermi line in the irreducible Brillouin zone. Numerical tests 
showed that well converged results can be obtained with about 300 Matsubara frequencies 
and 5 k-points along 1/8 of the Fermi line down to values for T c of ~ 0.002. Details of our 
method will be given elsewhere. 

The squares in Fig. 5 (joined by a solid line) show T c for r 3 -symmetry for doping 5 > 5 B o- 
The broken line exhibits T c if only the instantaneous part of the kernel is used. The dotted 
line corresponds to T c if the charge-charge contribution (last term in Eq.(l)) is dropped. 
Finally the dashed-dotted line describes T c if only the retarded kernel is taken into account. 

In order to understand the curves in Fig. 5 we consider the real part of the retarded 
kernel as a function of frequency, as shown in Fig. 4. Taking only the part between zero and 
uj 2 into account would yield rather large values for T c . To realize this one can decompose the 
real part into a large, constant repulsive part between uj = and uj = uj 2 plus the difference 
which is non-zero and attractive at low frequencies. Changing the cut-off from uj 2 to ou 1 using 
a pseudopotential description decreases strongly the effective repulsion yielding a large net 
attraction between and uj\ and thus a high value for T c . Model calculation show that the 
large negative part in the real part of ©3 above uj 2 is very harmful to superconductivity: 
Lowering first the cut-off from ~ 3\t\ to uj 2 increases the effective potential, reducing further 
the cut-off to uj\ decreases again the effective potential. The overall result is a rather 
modest attraction between uj = and uj — uj\. This explains the rather low values for 
T c calculated from the retarded part of alone, as shown by the dash-dotted line in Fig. 
5. An alternative consideration would start from the instantaneous part in leading to 
the transition temperatures shown by the broken line in Fig. 5. Using a realistic value 
for the energy unit of about 8000-K" (note that the usual effective hopping t exp corresponds 
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according to Eq.(l) to \t\/2) these values are already of the order of 100K and comparable 
to the experimental ones. Adding the retarded part of the kernel does not lead to a further 
increase in T c but rather to a slight decrease as illustrated by the squares in Fig. 5. Model 
calculations show that a purely attractive G rei would always increase T c . The observed 
lowering of T c must therefore be due to the repulsive part in G rei above U\. This part plays 
a role because the T c due to the instantaneous part alone is large enough to couple strongly 
the frequencies above and below ui\ in the gap equation. The doping dependence of T c is 
determined by various processes which compete with each other. For instance, the density of 
states decreases with doping which decreases the effective coupling. On the other hand , the 
energy scale for T c is in the case of the instantaneous term set by the effective band width 
which increases with doping. Fig. 5 shows that the net effect of these and other processes is 
to lower T c substantially with increasing doping for 5 > 5 bo- The eigenvalue A3 of the static 
kernel assumes according to Fig. 2 large negative values near 8 bo anticipating the incipient 
bond-order wave instability. The associated soft mode does not affect the instantaneous but 
the retarded part. Fig. 5 shows that the dashed line indeed increases steeply approaching 
5 bo from above. On the other hand little effects are seen in the squares representing the 
full calculation. The reason for this is that at the low T c 's of the dashed curve the relevant 
frequencies lie in the attractive region of Og et whereas for the high T c values of the squares 
these frequencies lie already in the repulsive part. From this one may conclude that neither 
soft modes nor large negative eigenvalues of the static kernel guarantee large T c values. In 
particular, there is no simple connection between T c and A3. 

The squares (connected by a dashed line) in Fig. 6 show T c as a function of the exchange 
constant J for a doping 5 = 0.17. T c approaches for J — > in agreement with previous 
findings^. Except at small values of J T c depends linearly on J. This is quite in contrast 
to the BCS-formula with an exponential dependence on J. The main reason for this is that 
the instantaneous term plays a major role and the momenta in it cannot confined to the 
Fermi surface in calculating T c . Eq.(53) is then a more appropriate formula for T c . It is, 
however, too simple to argue that 63 in that formula is rather independent on J and ~ 1/T 
for T > 0.01. The rather perfect linear behavior is the result of more subtle dependencies 
such as the J-dependence of the one-particle energies and the non-rigidity of the bands as a 
function of J. 

It has been argued0I3 that the nearest-neighbor Coulomb potential V is not neglegible 
in the cuprates. Adding this term 

H' = ^- £ XfXf, (63) 

ZJV <ij> 

p,q — l...N 

to the Hamiltonian Eq.(l) we have calculated T c as a function of V. The circles (connected 
by a solid line) in Fig. 6 present the result for J = 0.3 and 5 = 0.17. T c decreases strongly 
with increasing V and is extremely small for V > 3 J/4. This may be understood by looking 
just at the instantaneous term: Eq.(63) yields an additional term 2V(k — k') on the right- 
hand side of Eq.(43), canceling the J-terms exactly for V — J. T c would thus vanish if 
only the instantaneous term would be present. For V = J/2 H' cancels the charge term in 
Eq.(l). According to Fig. 6 this means a drop of T c by about a factor 5 compared to the 
value at V — which also agrees with Fig. 5. Such a big drop of T c due to the charge term 
of the t — J model seems very surprising because the latter can only produce effects ~ 5 2 
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and is therefore often omitted. Our calculation shows that this term cannot be neglected in 
a calculation of T c and increases T c substantially. 

Fig. 7 shows T c as a function of 5 for J = 0.3 and V = 0.15. Due to the Coulomb 
repulsion the T c values are rather low. The important frequencies in the solution of the gap 
equation are also low and mainly located in the attractive region of the retarded kernel, 
Adding the retarded to the instantaneous part thus increases T c . The two curves in Fig. 7 
demonstrate this effect. 

Figs. 8 and 9 show the influence of a second-nearest neighbor hopping term t' on T c . In 
these figures a lower value of 0.2 has been chosen for J similar as in RefB in order to have 
a reasonably low 5 bo an d a Van Hove singularity not too far away from optimal doping. 
Assuming a linear dependence of T c on J the the absolute values for T c are similar in the 
corresponding Figs. 5 and 8. T c is thus rather robust to changes in t'. T c decreases in Fig. 
8 somewhat slower than in Fig. 5 due to the larger density of states in the surroundings 
of the Van Hove singularity. The difference between squares and circles is also larger in 
Fig. 8 and decreases much slower with increasing doping. These effects are caused by the 
retarded part of the kernel which in Fig. 8 is less attractive below lu\ and more repulsive 
above uj\ as shown in Fig. 2 of Ref.i (Note that the curve in that Figure includes a constant 
instantaneous contribution of -0.29). The Van Hove singularity is located near 5 = 0.28. 
There is nearly no effect of the van Hove singularity on T c for the same reasons as in the 
above case of the incipient bond-order wave instability. (The curves in Fig. 3 of Ref.i were 
calculated without the charge-charge term and, erroneously, without the factor 1/2 in the 
instantaneous part. Correcting this error amounts essentially to lower the curves in that 
Figure by about a factor 4). Fig. 9 should be compared with Fig. 6. T c decreases in both 
cases strongly with increasing V and practically vanishes for V > J. 

We have also searched for superconducting instabilities with order parameters which are 
odd in frequency. The sum over Matsubara frequencies is always zero in this case so that 
the constraint of having no double occupancies of sites at the same time is automatically 
fulfilled for the two particles of the Cooper pair. There exists no static approximation for the 
kernel in this case. Possible instabilities are again determined by the linearized gap equation 
Eq.(40) where S an and have to be projected on the odd frequency parts. The transition 
temperature is determined by the condition that the determinant of a matrix consisting 
essentially of the kernel and the unity matrix is zero. In Fig. 10 we have plotted this 
determinant as function of the temperature for each of the five irreducible representations 
using J = 0.3, t' = 0., and 5 = 0.17. The Figure clearly shows that none of the curves 
tends to zero in the investigated temperature interval ruling out any odd frequency pairing 
instability with a T c larger than ~ 0.002. This also can be seen directly from Fig. 11 
where the real part of Q3(u> + if]) is shown for odd frequency pairing with Ta-symmetry 
corresponding also to triplet pairing. The effective interaction is repulsive up to energies ~ t 
ruling out an instability towards superconductivity in this channel. 

D. Conclusions 

Sections II and III demonstrate the feasibility of developing a perturbation expansion for 
the t — J model in terms of X-operators and obtaining explicit expressions for the leading 
contributions to the anomalous self-energy of physical electrons using an 1/N expansion. The 
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nonapplicability of Wick's theorem to X-operators does not cause any problem in such an 
approach. One also should note that the more involved and sophisticated parts of sections II 
and III deal with contributions which are far beyond those considered in previous treatments. 
For instance, the slave boson 1/N result! for the anomalous self-energy corresponds in our 
approach just to the inhomogenous term in the integral equation Eq.(32_). Similarly, only 
part of the first contribution to the kernel, O^, was calculated in Ref.Q using a diagram 
technique for X-operators. This shows in our opinion that the employed functional approach 
is at least as suitable as other approaches to treat highly correlated Fermi systems. Taking 
also the numerical results of section 3 into account our main conclusions can be summarized 
as follows: 

a) Our explicit expression for the 0(1/N) anomalous self-energy is clearly different from 
the corresponding expression of the slave boson theory. In particular, the largest contribution 
to the retarded kernel comes in the present approach from the anomalous part of the vertex 
and has no analogue in the slave boson approach. The presented expressions show for the first 
time in an explicit way that the 1/N expansions are really different in the two approaches. 
This difference can be traced back to different Hilbert spaces and different enforcement of 
the constraint. Even the order parameters for superconductivity are not related in a simple 
way: The leading slave boson order parameter involves necessarily (small) violations of local 
constraints in order to be nonzero whereas such violations are ruled out in our approach. 

b) The kernel of the linearized gap equation consists of an instantaneous and an re- 
tarded part and both are similar in magnitude at low frequencies and for momenta near 
the Fermi surface. We found that there are superconducting instabilities in each symmetry 
channel and for all dopings. The true ground state thus never describes a Fermi liquid but 
a superconductor similar as in the weak coupling casell. However, these instabilities are 
in general extremely weak leading to academically low transition temperatures. The only 
clear and robust exception is the d-wavelike symmetry where a strong instability towards 
superconductivity occurs. Odd symmetry pairing mechanisms turned out to be very weak 
and can be ruled out as a mechanism for high-T c superconductivity in our model. 

c) In the case of T 3 symmetry the real part of the retarded kernel is weakly attractive at 
low frequencies on an energy scale of J or a fraction thereof and strongly repulsive at higher 
frequencies whereas the instantaneous part is attractive. Solving numerically the linearized 
gap equation the obtained transition temperatures T c are for N = 2 of the order of 0.01|t| and 
thus in principle large enough to account for the phenomena of high-T c superconductivity. 
It is interesting to note that the Hubbard model at small or intermediate couplings also 
shows d-wave superconductivity with similar values for T c .il The instantaneous term is 
instrumental in getting these large values for T c : First, its cutoff is given by the effective 
band width and thus in general larger than J. Secondly, due to the large T c , the solution 
of the gap equation involves large frequencies where the retarded term is strongly repulsive. 
As a result the retarded term is of less importance because attractive and repulsive parts 
of it cancel each other to a large extent. The dominance of the instantaneous part and 
the presence of two cutoffs lead to strong deviations from BCS-behavior. For instance, T c 
depends linearly and not exponentially on J except at very small values for J. 

d) T c is rather insensitive to the addition of a second-nearest neighbor hopping term t' 
and to a Van Hove singularity. The latter can be understood by noting that the solution 
of the gap equation involves momentum and frequency averages in the instantaneous and 
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retarded part, respectively, so that singularities in the density of states are washed out. T c 
depends, however, sensitively on a nearest-neighbor Coulomb repulsion V and becomes very 
small if V is substantially larger than J. 

e) Our calculations are limited to dopings larger than 5bo(~ 0.14 for J = 0.3) where 
an instability towards an incommensurate bond-order wave of d-symmetry occurs. The as- 
sociated soft mode causes A3 — > —00 for 5 — > 5bo whereas T c is nearly unaffected. In the 
underdoped regime 5 < 5bo we expect a competition between bond-order and antiferromag- 
netic fluctuations which is beyond the leading order of the 1 /N expansion considered in this 
investigation. 
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FIGURES 



FIG. 1. Dependence of various contributions to the total kernel on the second momentum 
k' along the Fermi line for a fixed first momentum k = (2.465,0.309); a) contributions Q^\@^; 
b) contributions 9 (3) , 9 (4) , and the total 9. 

FIG. 2. Lowest eigenvalues A« of the static kernel 9 for the five representations Tj of C± v as a 
function of the doping 5. 

FIG. 3. Negative imaginary part of the d-wave projected kernel 93 {u + irf) as a function of 
the frequency ui using r\ = 0.005. 

FIG. 4. Real part of the d-wave projected kernel 93(0; + irj) as a function of the frequency ui 
for 77 = 0.005. Inset: The same for small frequencies using i] = 0.002. 

FIG. 5. Transition temperature T c (with \t\ as energy unit) for T3 pairing using the total 
kernel (squares), the instantaneous part (circles), the retarded part (dimaonds) of the kernel, and 
the total kernel without charge-charge term (triangles). 

FIG. 6. Transition temperature T c for T3 pairing versus V (circles) and versus J (squares). 

FIG. 7. Transition temperature T c for T3 pairing for V = 0.15 using the total kernel (squares) 
and the instantaneous part of the kernel (circles). 

FIG. 8. Transition temperature T c for T3 pairing for J = 0.2 and t' = —0.35 using the total 
kernel (squares) and the instantaneous part of the kernel (circles) . 

FIG. 9. Transition temperature T c for J = 0.2, t' = —0.35, 5 = 0.23 as a function of V using 
the total kernel (squares) and the instantaneous part of the kernel (circles) . 

FIG. 10. Determinant associated with the gap equation for odd frequency pairing with sym- 
metry Ti as a function of temperature T. 

FIG. 11. Real part of the odd frequency, T3 kernel 93(0; + irf) as a function of the frequency 
uj for i] = 0.005. 
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